CHINESE JOURNAL 
OF ENGINEERING MATHEMATICS 


June 2010 Vol. 27 No. 3 


Article ID:1005-3085(2010)03-0534-15 


Mathematical Modeling and Dynamics of HIV 
Progression and Treatment* 


WANG Yan, ZHOU Yi-cang 
(Department of Applied Mathematics, Xi’an Jiaotong University, Xi’an 710049) 


Abstract: This article presents a brief review on the modeling and dynamics of HIV infection in 
vivo. We introduce typical mathematical models on the interaction between CD4* T- 
cells and virus particles, as well as the drug therapy. We focus on theoretical results 
and simulations of ODEs, DDEs, integro-differential equations and impulsive differential 
equations. Parameter values of those models are collected in a table for reader’s reference. 
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1 Introduction 


AIDS (Acquired Immunodeficiency Syndrome) is one of the most serious, deadly diseases 
in human history. AIDS, caused by the human immunodeficiency virus (HIV), has spread in 
successive waves in various regions around the globe. The estimated number of persons living 
with HIV worldwide in 2008 was 33.4 million. The number of newly infected with HIV in 2008 
is 2.7 million”. 

HIV needs bonding to the specific host cell receptor. The receptor is a CD4 antigen in the 
membrane of T helpers lymphocytes. When the cell is activated by an antigenic impulse, viral 
protease finishes the process of releasing the virus. These new virions immediately attack next 
CD4* T-cells and leads to their destruction. The immune system is exhausted by acceleration 
of the CD4+ production and destruction, and soon the destruction predominates. The deep 
CD4t count decreasing (under 200/mm?) is strongly associated with the risk of development 
of opportunistic infections. 

Although treatments for AIDS and HIV can slow AIDS progression, there is currently no 
vaccine or cure. The fact that HIV replicates rapidly, producing on average 10?° viral particles 
per day, leads to the realization that HIV evolves so rapidly that treatment with a single 
drug is bound to faill?l. Highly Active Antiretroviral Therapy (HAART) is the recommended 
treatment for HIV infection. HAART combines three or more anti-HIV medications in a daily 
regimen. The drugs are unable to eradicate virus from the body, but they can minimize viral 
replication, restore immunological function, improve quality of life and reduce morbidity and 
mortality. These drug cocktails generally consist of reverse transcriptase (RTIs) and protease 
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inhibitors (PIs). RTIs can block the infection of target T-cells by infectious virus, and PIs can 
prevent HIV protease from cleaving HIV polyprotein into functional units, causing infected cells 
to produce virus particles that are noninfectious. There are another two new agents recently: 
entry inhibitors and integrase inhibitors. 


Viral mathematical models can provide insights into the dynamics of viral load in vivo, 
and may play a significant role in the development of a better understanding of diseases and 
various drug therapy strategies against them. Perelson and Nelson!*), Nowak and May!) provide 
excellent reviews and many more citations on HIV treatment. This review paper will give a 
brief introduction to modeling and dynamics on the interaction between HIV and CD4t T-cells. 
The main attention is paid to mathematical analysis of interaction between T-cells, antiviral 
immune system responses, drug therapy and viral load. Furthermore, theoretical results and 
methods of viral microcosmic models are included in paper'4!. For convenience, we use T-cells 
instead of CD4+ T-cells in the rest parts of the paper. 


2 Basic models with logistic T-cells growth and immune response 


2.1 The basic model and the extension 


We begin with a very simple model, which captures essential features of HIV. The basic 
model of viral dynamics has three variables: uninfected T-cells concentration at time t, T(t); 
infected T-cells concentration at time t, T*(t); free viral particles concentration at time t, V(t). 
Uninfected T-cells are produced at a constant rate, A, and die at a rate dT’. Free virus infects 
uninfected T-cells to produce infected T-cells at rate kTV. Infected cells die at rate 6T*. New 
virus is produced from infected cells at rate NôT* and dies at rate cV. The average life-times of 
uninfected T-cells, infected T-cells, and free virus are given by 1/d, 1/6, and 1/c, respectively. 
The average number of virus particles produced over the lifetime of a single infected T-cell (the 
burst size) is given by N. These assumptions lead to the following differential equations 


T (t) = à -dT —kVT, 
T* (t) = kVT — 6T", (1) 
V(t) = NOT* — cV. 


Model (1) has been studied in [3,5-7]. (1) has an uninfected steady state Eo(A/d,0,0). The 
crucial quantity is the basic reproductive ratio, Ro = (kKNA)/(cd). Ro is the average number of 
newly infected cells that arise from one infected cell when almost all cells are uninfected. When 
Ro < 1, the uninfected steady state Ep is locally asymptotically stable; when Ro > 1, Eo is 
unstable while the infected steady state E = (A/(dRo), cd(Ro — 1)/(k Nô), d(Ro — 1)/k) exists 
and is locally asymptotically stable. 

Considering the density-dependent proliferation rate of T-cells!®l, logistic growth has been 
included in other models to describe the growth of healthy T-cells!#*14]. In 2003, De Leenheer 
and Smith!®! summarized a competitive model reflecting the interaction between T-cells and 
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virus particles, given by following equations 


T(t) = f(T) — kVT, 
T* (t) =kVT — ôT*, (2) 
V'(t) = NôT* — cV —ikVT, 


where 
T : ; 
f(T) =\-dT +rT(1--—), or f(T)=A-—dT, and i=0, or i=l. 


r is the growth rate (r > d in general), and Tmax is the carrying capacity. Using the theory of 
competitive (and cooperate) systems by Hirsch and Smith (see [9] and the references therein), 
those authors obtained the basic reproduction number Ro = KTy(N —i)/c (To is the normal 
target cell level at the uninfected steady state) for model (2) and got following results. 

1) If Ro <1, the disease-free steady state is globally stable and the virus is cleared, which 
does not depend on the form of f(T). 

2) IfRg >1and f(T) = A—dT, then a chronic disease steady state exists which is globally 
asymptotically stable under certain conditions. 

3) If Ro > 1 and 


f(T)=à-dT+rT(1- r ), 


there exists an orbitally asymptotically stable periodic orbit, which attracts almost all solutions 
under suitable conditions. 

It was shown tnhat particular choices for f(T) may lead to different qualitative behavior. 
For example, for f(T) = A — dT, the chronic disease steady state, if it exists, is always locally 


asymptotically stable, while for 


f(T) =A-dT+rT(1- K | 


this steady state may be unstable and sustained oscillations may occur'!. This sensitivity of 
the behavior to f(T), in particular, calls for a better understanding of the mechanisms of T-cells 
proliferation. 

Wang and Ellermeyer studied a model with a full logistic proliferation term rT(1 — wr) 
and got the stability of infected steady state depending on the size of the T-cells proliferation 
rate, Culshaw and Ruan chose N as a bifurcation parameter and obtained a restriction 
Nerit on the number of viral particles released per infectious cell in order for infection to be 
sustained!!4], Wang and Lil? obtained the global stability of the positive equilibrium by the 
compound matrices method. 

2.2 Immune response 

For most virus infections, cytotoxic T lymphocytes (CTLs) play a critical role in antiviral 
defense by attacking virus-infected cells. It is believed that they are the main host immune 
factor tnhat limits the extent of virus replication in vivo and thus determines virus load!®!. The 
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mathematical model with CTLs is a system of differential equations 


T (t) =\—dT —kVT, 


T* (t) = kVT — 6T* — qT*C, T 
3 
V(t) = N6T* — cV, 


C(t) = g(T,T*,C) — dcC, 


where C(t) denotes the concentration of CTLs at time t, and dg is the death rate of CTLs. 
Infected T-cells die at a per capita rate 6 and are killed by CTLs with a rate gT*C. The 
function g(T,T*,C) describes the rate of activation of the immune response (CTLs). It has 
been suggested that g(T,T*, C) can have one of the following three forms!® 913] 


g(T,T*,C) =pT*, g(T,T*,C)=pT*C, g(T,T*,C)=pTT*C, 


where p is the rate of stimulation of CTLs. It is still unknown which g(T,T*,C) is better 
to describe the activation of CTLs. Based on this simple model and simulations, Nowak and 
Bangham explained that CTL responsiveness determines virus load, but there may be no obvi- 
ous correlation between virus load and the abundance of antiviral CTLs!®]. They showed that a 
better indicator of CTL responsiveness is the equilibrium virus load, rather than the abundance 
of virus-specific CTLs. 


3 Time delay 


It is well known that delay differential equations (DDEs) exhibit much more complicated 
dynamics than ODEs since a delay may cause the stability change and oscillation. Actually, time 
delay occurs in the progression of HIV infection. In studying the viral clearance rates, Perelson 
et all assumed that there are two types of delays, pharmacological delay and intracellular 
delay, which occur between the administration of drug and the observed decline in viral load. 
The pharmacological delay is the time from drug administration to the entry of the virus within 
the target cell. The intracellular delay is the time required for an infected cell to replicate (the 
time interval that the virus binds to the receptors of the target cell to the time the first virion 
is produced from the same target cell). 

3.1 Discrete delay 

In 1996, Herz et all] used a discrete delay to model the intracellular delay in the basic HIV 
model and showed that the incorporation of a delay would substantially shorten the estimated 
half-life of the free virus. Nelson et alë] found that the predicted rate of decline in plasma virus 
concentration depends on three factors: the death rate of virus producing cells, the efficacy of 
therapy, and the length of the delay if an intracellular delay and drug therapy are included. 

Culshaw and Ruan extended the basic model to incorporate the logistic growth and an 
intracellular delay!!]. They demonstrated that the stability of the infected steady state is not 
affected by the modification with realistic parameter values. However, they did find conditions 
that lead to a Hopf bifurcation if parameter intervals are beyond certain intervals. We will 
show a delay model and its dynamics in Section 4. 
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3.2 Continuous delay 

We often use gamma distribution to represent multistage processes. We can assume that 
the probability distribution of infected T-cells death is given by a gamma distribution, which 
suggests that the death of infected T-cells only occurs after a number of subprocesses are 
completed], A gamma distribution function is used to describe a continuous delay between 
infection and viral production[16, 17]. Nelson and Perelson!'®! constructed an integro-differential 
equation and showed that the value for ô estimated from patient data was directly influenced 
by the assumed variance and mean of the gamma distribution, and concluded by presenting a 
general result on the stability of a set of DDEs. Mittler et all!”] showed that good estimates for 
free viral clearance rates, infected cell death rates, and parameters characterizing the gamma 
distribution can be obtained. They demonstrated that it is possible to incorporate distributed 
intracellular delays into existing models for HIV dynamics and to use these refined models to 
estimate the half-life of free virus from data. 


4 Drug therapy 


4.1 Combination therapy 

Research has shown that using combinations of drugs is a better treatment strategy than 
using only one or two drugs!?]. We focus on HIV models with the combination of RTIs and 
PIs in order to reduce the amount of virus in their bodies. Let ny and Tip be the drug efficacy 
of RTIs and PIs, respectively. We assume the drug efficacy is less than 100% in the following 
models. 

Based on system (2), De Leenheer and Smith!) studied the following model with constant 
drug therapy 

T(t) = f(T) — k(1 — nye) WiT, 


T* (t) =k(1—n»)V/T — 6T*, 
V; (t) = (1 — np) NOT* — ikV;T — cV}, 
Vy (t) = NpNOT* — cVyr. 


The basic reproduction number of (4) is R§ = kTp[N(1—nre)(1—np) —i]/c. Similar conclusions 
as those for system (2) were obtained. We continue to investigate the interaction between T-cells 
and virus load with drug therapy. We formulate the model with HIV pathogenesis within-host, 
adding the combined effects of drug therapy (RTIs and PIs), an intracellular delay and logistic 
growth!4] 


T(t) =—dT(t)+ rr(e)(A To ) ky(1 — na) V (ÐT (E), 
T* (t) = ka(1 — nr) Vi(t — T)T(t — 7) — 6T*(t), 

V; (t) = (1 — np) NST” (t) — cVr(t), 

Vyr(t) = np NST” (t) — Vyr (t). 


(5) 


The parameter k; represents the rate of infection of T-cells with free virus, and k4 is the rate 
at which infected cells become actively infected. We assume that dTmax > A, the number of 
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T-cells will decrease if it surpasses Tmax. The dynamics of system (5) is fully determined by 
the first three equations since 


t 
Vynr(t) = e= ( Vy 1(0) + | e° nyN6T"(r)dr ). 
0 
The initial conditions are 
T(0)=91(0)>0, T*(0)=y2>0, Vi(@)= (0) 20 for 6 € {—7,0], (6) 


where ¥ is a given constant, and y1, Y € C({-7,0], R+) with Ry = (0, +00]. It can be verified 
that if (T(t), T*(t), Vr(t)) is a solution with initial conditions (6), then T(t) = 0, T*(¢) > 
0, Vr(t) > 0 for all t > 0. Let 


c 
Ne = 1 — (1 — Nnrt)(1— np) Nerit:=1— ANI 


Then Ne < Nerit ensures the existence of the infected (positive) steady state E. 

We obtained the transcendental characteristic equation F} (€) of the linearized system at 
the infected steady state Æ. When 7 = 0, we got the asymptotic stability condition of E by 
Routh-Hurwitz criterion. When 7 > 0, using Kuang’s theory of DDEs!!8], we proved that F, (£) 
can not have a purely imaginary root under certain conditions, and the infected steady state 
E is locally asymptotical stable. Under other conditions we concluded that the characteristic 
equation F; (£) has a pair of purely imaginary roots tiwp. We proved the occurrence of a Hopf 
bifurcation at the positive equilibrium E. Let €(7) = x(r) + iy(r) be the root of Fi(€), and 


S= sign{ $ (Ree()} Ea sign{Re( T yes 


Differentiate F\(€) with respect to 7 and combine the above formula, we can obtain sufficient 


conditions to satisfy S > 0 (the transversality condition). 

By the theory of DDEs, sufficient conditions are established for the local asymptotic stability 
of the uninfected steady state and the infected steady state. Furthermore, by constructing 
Lyapunov function, we obtained the global asymptotic stability of uninfected steady state when 
Ne > Nerit- The influence of the time delay on the stability of equilibrium states is discussed. We 
showed that the local stability of the uninfected steady state is independent of the size of the 
delay; on the other hand, we proved that increasing the delay can destabilize the infected steady 
state leading to a Hopf bifurcation of periodic solutions. This clearly shows the importance of 
time delay on HIV dynamics under treatment. 

In simulation, we have parameters À = 10 cells mm~*day~+, ô = 0.26 day~* and Tmax = 
1500 mm~® fixed!4], and choose three sets of data: nre = 0.4, np = 0.3; mrt = 0.5, np = 0.55 
and nre = 0.6, np = 0.8; then ne = 0.5800, ne = 0.7750 and Ne = 0.9200, respectively. The 
numerical simulations show that increasing the combination drug efficacy can increase the total 
number of healthy and infected T-cells and decrease the number of the infectious virus (see 
Figure 1). In other case, we choose n,e = 0.6, np = 0.7, then ne = 0.8800 and Ne < Merit: 
Comparing to Culshaw and Ruan!!4), we found a Hopf bifurcation in realistic parameter values 
(see Figure 2). Similar conclusion can be obtained if we use the full logistic growth of healthy 
T-cells in model (5). Our model is an extension of Culshaw and Ruan!” 
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Figure 1: Local asymptotic stability of the infected steady state E. d = 0.007 day™*, r = 0.03 day, 
kı =2.4 x 1075 mm3day~!, ka = 2.0 x 107° mmday~!, N = 2500, c= 2.4 day™!, T = 1.5 

days, and nc changes from 0.5800 to 0.9200. As ņ« increases, the total number of the healthy 

and infected T-cells (a) increases dramatically, while the number of the infectious virions (b) 


decreases substantially 
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Figure 2: Periodic solutions bifurcated from the infected steady state Ē. d = 0.03 day™!, r = 0.95 
day~+, kr = 0.0027 virions mm?day™*, ka = 0.0020 virions mmĉday™', N = 50, c=3 


day! and 


T = 1.5 > To & 1.3395 


The constant efficacy of RTIs and PIs is only a simple assumption in modeling HIV treat- 
ment. The efficacy of antiviral treatment varies with time and other factors. Mathematical 
models have been developed to study the pharmacokinetics of drug therapies!!9.20l, Huang et 
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all?°] formulated a model to illustrate the relationship between antiviral response and phar- 
macokinetics, time-varying adherence and drug resistance by simulation experiments. They 
discussed the properties of the viral dynamics and studied how time-varying treatment effica- 
cies affect antiviral responses, specifically viral load and T-cell count. 

In the long term of HIV treatment, the drug concentration varies during the dose intervals. 
Once a dose is administered, drug concentration gets high and reaches a peak value. As time 
passes by, the concentration gets lower. When another dose is administered, drug concentration 
repeats the same cycle. It is natural to assume that drug concentration varies periodically, which 
results in periodic variation of the drug efficacy. Yang and Xiao!) formulated following simple 
model with periodic drug efficacy and investigated the interaction dynamics of T-cells and virus. 


T (t) =A — dT — k(1 — na (t) VT, 
T* (t) = k(1 — na lt)) VT — ôT*, (7) 
V (t) = N6T* — cV. 
Yang and Xiao got a threshold value of (7) between the extinction and the uniform persis- 
tence of the disease by applying the persistence theory. The disease goes to extinction if the 
threshold value is less than unity, while the disease persists if the threshold value is larger than 
unity. There exists a positive periodic solution which is globally asymptotically stable. Their 
conclusion extended the classic results for the basic autonomous model. 
Based on Dixit and Perelson!!9] model, we extend system (5) to include pharmacokinetics 
of drug therapies. 
T (t) =A— dT (t) — (1 — nalt) kV OTE), 
T* (t) = (1 — nrelt — 7))kV7(t ~ T)T(t — Te — ôT* (t), 
Vi (t) = (1 — np(t)) NST* (t) — cVi(t), 
Vy 7(t) = np(t)N6T*(t) — cVivr(t). 


(8) 


The time-varying parameters n,,(t) and n,(t) are the RTJs and Pls drug efficacies at time t, 
respectively. The challengeable problem of (8) is to determine the expressions of n,,(t) and 
Np(t) since there is no good way to get them by laboratory experiment. We are trying to 
collect parameter values and treatment data to investigate the dynamics of (8) by numerical 
simulations. 

4.2 Mutation and drug resistance 

HIV does not always make perfect copies of itself. With billions of virions being made 
every day, small and random differences can happen. The differences are called mutations. 
Mutations can keep the drugs from working, leading to the reduction in effectiveness. This 
is called drug resistance. The drug resistance is closely related to the mutation and poor 
adherence. Adherence means taking your medications on time, taking the prescribed dose, and 
taking them the correct way (with or without food, for example). 

The emergence of drug resistance is one of the most prevalent reasons for treatment failure in 
HIV therapy. There are a lot of mathematical models concerning the population dynamics, HIV 
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treatment and drug resistance!®.22:23], Recently, Rong et all2?] assumed that the drug-sensitive 
and resistant strains differ by a single mutation. They first showed that for chronically infected 
HIV patients the drug resistant strain exists before the initiation of antiretroviral therapy 
(ART). The model is the following system of differential equations 


T (t) =A — dT — ks V,T — k, VT, 

T; (t) = (1 — u)ks VT — OTs, 

V; (t) = N,OTs — CVa, (9) 
T(t) = Uks VT + k, V,T — ôT p, 

V, (t) = NOT, — eV, 


where, T;(t) and T, (t) are the concentrations of T-cells productively infected by drug sensitive 
virus and drug-resistant virus, respectively; V(t) and V,(t) represent the concentrations of 
drug sensitive and drug-resistant virus, respectively. k, and k, represent the rate constants 
at which uninfected cells are infected by drug sensitive and drug-resistant virus, respectively. 
u(0 < u < 1) is the rate at which cells infected by the drug sensitive virus mutate and become 
drug-resistant during the process of reverse transcription of viral RNA into proviral DNA. Both 
types of infected cells are assumed to have the same death rate 6. They supposed that the drug 
sensitive and resistant strains differ in their burst sizes, N, and N,, while they have the same 
virion clearance rate c. It should be noted that the backward mutation from drug-resistant to 
drug-sensitive strain is neglected since the wild-type virus dominates the population before the 
initiation of therapy!”]. They assumed both the infection rate and burst size of resistant strain 
are less than those of wild-type strain, i.e., kp < ks and N, < Ng. 

Those authors aimed to study the virus dynamics as well as the development of drug- 
resistant strains. They gave the basic reproductive numbers of the wild-type strain R, and the 
drug-resistant strain R,, respectively, Rs = (ksNsX)/(dc) and R, = (k,N,A)/(dc). They got 
following results: 

1) The infection-free steady state Eo is locally asymptotically stable if Ra < 1/(1 — u} and 
Rr < 1, and it is unstable if Ra > 1/(1 — u) or R, > 1. 

2) The steady state with only drug-resistant virus, Æp, exists if and only if Ry > 1. It is 
locally asymptotically stable if R, > (1 — u)/Rs and unstable if R, < (1 — u)/Rs. 

3) The coexistence steady state E, exists and is locally asymptotically stable if and only 
if Rs > 1/(1 — u) and R, < (1 ~ u)/Rs. 

With parameter values \ = 104 ml`!day™!, d = 0.0lday~!, k, = 2.4 x 1078 ml day™}t, 
kr = 2.0 x 1078 ml day~1, u = 3 x 1075, 6 = iday~!, Ns = 3000, N, = 2000 and c = 23day—!, 
Figure 3 shows simulation results for the T-cells count and viral loads of both the wild-type and 
drug-resistant strains before treatment. The wild-type virus dominates the virus population 
before the initiation of ART. 

4.3 Impulsive drug effects 

Recently, Smith and Wah!!?4,25] described the drug concentrations during HIV therapy using 
impulsive differential equations, considering the dynamics of T-cells interacting with free virions, 
reverse transcriptase inhibiting drugs and protease inhibiting drugs. In [24], they divided T- 
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cells into six classes and found that insufficient dosing of either drug corresponds to high viral 
load and a large population of infectious T cells. They further predicted that, in the absence 
of physiological limits on tolerable drug concentrations, sufficiently frequent dosing with the 
RTI alone could theoretically maintain the T-cell count close to the uninfected case. However, 
for frequent dosing of the PI alone, the limiting T-cell populations may not be enough to 
maintain the immunity. Furthermore, frequent dosing of both drugs has the same effect on the 
T-cell population as frequent dosing of the single RTI. Those two drugs can have fundamentally 
different effects on the long-term dynamics and the RTI, in particular, plays a crucial role in 
maintaining the immunity. In [25], different regimes were classified according to whether the 
drug efficacy is negligible, intermediate or high. Two strains of virus were considered: a wild- 
type strain that can be controlled by both intermediate and high drug concentrations, and a 
mutant strain that can only be controlled by high drug concentrations. It was discussed if 
dosing schedules and concentrations of preventative drugs facilitate or prevent the emergence 
of drug resistance. The estimates of a range of dosages and dosing schedules which would, if 
physiologically tolerable, theoretically eliminate free virus in this system were provided. Their 
results predicted that decreasing the interval between doses was more effective than increasing 
the dose to contro} viral load. 
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Figure 3: Simulation of uninfected T-cells and viral loads of both wild-type and drug-resistant 
strains for the pretreatment model (9). The basic reproductive numbers for two strains 
are R, = 3.13 and R, = 1.74. Before treatment, both strains of virus coexist. However, 
the drug-resistant strain remains at a very low level after an initial increase and the 


wild-type virus dominates the population 
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5 Conclusion and discussion 


We have presented a number of dynamical models and results on the study of HIV dynam- 
ics. We have tried to demonstrate how mathematical modeling can help us to understand HIV 
pathogenesis. The interaction of T-cells and virus, time delay, drug therapy and drug resis- 
tance are discussed in different models. Mathematical modeling and computer simulations are 
powerful tool to improve the design of AIDS clinical trials, to study the pathogenesis of HIV 
infection, to estimate the parameter values, and to evaluate new treatment strategies!” 

Most of the mathematical models can be extended to a more general system 


è = f(x)(or f(z,y))-—2z, y=az(ora(t—r)2(t—7))-y, ż=y-z, 


which, for example, can be used describe the interaction of caterpillars x, parasitic wasps z and 
parasitized caterpillars y. Those models and results can be applied in other settings beyond 
the HIV dynamics. It is hopeful to obtain more mathematical results and richer dynamics for 
the model with a general function f. 

Although there are plenty of mathematical models and simulations, there still a long way 
for us to understand the HIV pathogenesis completely. There are many interesting and chal- 
lengeable topics in modeling and dynamics of HIV infection and treatment. 

5.1 Age and infection age structures 

Age and infection age are important factors in HIV infection and disease progression. Age- 
structured models have been widely used to study the epidemiology of infectious diseases. 
Nelson et all?! developed an age-structured model of HIV infection and showed that the age- 
structured model is a generalization of ODE and DDE models. Age-structured models have 
greater flexibility that may help us to understand the underlying mechanism of HIV infection 
and treatment better?27]. Infection age-structured model may be more appropriate in describing 
the disease progression and treatment since HIV takes about 10 years, on average, to progress 
from initial infection to AIDS, T-cell count and viral load varying a lot in different stages of 
infection. Infection age-structured model may give better simulation and assessment result of 
a treatment strategy. 

5.2 Intermittent treatment 

Kirschner and Webb!23] showed that a possible strategy, intermittent treatment, may have 
advantage in reducing side effects or resistance. Intermittent strategy is to treat with chemother- 
apy followed by interruptions in the treatment during which either a different drug or no 
treatment is administered. Recently, the Strategies for Management of Antiretroviral Ther- 
apy (SMART) Study Group!*! found that the intermittent ART guided by the CD4t count 
significantly increased the risk of opportunistic disease or death from any cause, as compared 
with continuous ART, largely as a consequence of lowering the CD4* cell count and increasing 
the viral load. Intermittent ART does not reduce the risk of adverse events that have been as- 
sociated with ART. However, because of the poor adherence in reality, or by other interruption 
of treatment, it is necessary to investigate the intermittent treatment. 

5.3 HIV co-infection with other diseases 

Kirschner!2°] indicated that co-infection may indeed play a dramatic role in diseases. He 
found that the T-cell populations are lower in the presence of both M. tuberculosis and HIV 
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than in the case of infection with HIV alone. Also, the viral load and M. tuberculosis population 
are higher in the co-infected patient, than the single-pathogen infection cases. 

Co-infection with hepatitis C virus (HCV) and HIV is common in certain populations. 
Among HCV(+) persons, 10% are also HIV(+), and among HIV(+) persons, 25% are also 
HCV(+). Many studies have shown that in intravenous drug users, co-infection prevalence 
can be as high as 90-95%. There is increasing evidence supporting the concept that people 
infected with HIV have a much more rapid course of their hepatitis C infection®%. HIV can 
also co-infection with Hepatitis B virus (HBV). 

5.4 Latent infection 

Although the HAART is extremely effective in reducing the viral load, which resulted in 
a substantia] reduction of morbidity and mortality, an HIV infected individual has not got 
complete recovery by the prolonged HAART. HIV can have latent infection in resting memory 
T-cells, which forms a cellular reservoir. The nature of memory T-cells have them remain in 
the resting state in the presence of HAART for a long time. The viral may rebound when 
HAART is withdrawn since reservoir releases new virus. Rong and Perelson presented a review 
on viral persistence, the latent reservoir, and mathematical models developed to explore their 
relationships!®2!. Interested readers can have more information from their paper. Banks et al 
formulated a model to describe the pathogenesis of HIV infection including certain important 
features, such as the reservoir of latently infected cells!33]. That more biological models may 
admit multiple stable off-treatment equilibria, and can exhibit the phenomenon of transient 
viremia. They showed that the model provides reasonable fits to the patient data and, moreover, 
that it exhibits good predictive capability. They demonstrated that parameter values obtained 
for most clinical patients do not admit multiple stable off-treatment equilibria. 

5.5 Parameter estimate and model fitting with data 

In the development of mathematical models for the progression of HIV infection and treat- 
ment, the model dynamics should be in agreement with the mechanism and key features of 
HIV infection. In application of those models to HIV treatment the parameter values should 
be determined, and the model prediction should match with data. Adams et al discussed 
mathematical and statistical ideas relevant to structured treatment interruptions, including 
parameter estimation!34!. A nonlinear model was used to fit clinical data for HIV patients 
undergoing treatment interruptions. After the parameter estimation by a statistically-based 
censored data combined with inverse problem techniques the comparison between the model 
simulation result and observation data was demonstrated by many figures!?5). 

Although mathematical models and theoretical results have done great help to understand 
HIV infection and treatment, the model assumption and parameter estimation method should 
be improved to make the model and the application more biological and realistic. If we use 
model with constant parameters to describe the long term dynamics of HIV infection and 
treatment, it is usually found that the T-cell count and the viral load will tend to equilibrium 
states. The phenomenon tending to equilibrium states of those models does not match the 
biological observation. The T-cell count (the viral load) of an HIV infected individual will be 
very low (very high) after becoming AIDS without treatment. We are working to develop a 
model with time dependent parameters to simulate the complete progression of an HIV infected 


546 CHINESE JOURNAL OF ENGINEERING MATHEMATICS VOL. 27 


individual. We determine those time-dependent parameters by comparing the model simulation 
result with the well known curves of T-cell count and viral load. After those parameters are 
confirmed, we can use the model to give more biological prediction on the disease progression 
and treatment. 

There are also other challenging topics in HIV infection and treatment. 1) HIV treatment 
in a specific country: in a specific country, the environment, HIV transmission, the regimes of 
treatment, and the combination of ART may be different from other countries; 2) Time depen- 
dent drug efficacy: it is more realistic assumption for the interaction between pharmacokinetics 
and pharmacodynamics. We may obtain the expression among the regime treatment, drug 
efficacy with plasma concentration by experimentation; 3) Global dynamics: We do not have 
effective method to obtain satisfactory results on global dynamics. For most of the models, 
including our modell], the global stability is still an open question. New theory and method 
need to be established. 


6 Appendix 


We have searched many references to determine the range of parameter values. We list them 
in the appendix so that our readers can consult easily, see Table 1 and all the parameters cited 


from!“ 
Table 1: List of parameters 
Parameter Definition Range of Parameters 
T-cells source term 0-10 mm~*day~! 
Death rate of healthy T-cells 0.007-0.1 day~? 
r Growth rate of T-cells 0.03-3day~* 
ky Viral infectivity rate 0.00025-0.5 mm*day~* 
ka Rate infected cells becomes active ka/kr %1 
ô Death rate of infected T-cells 0.2-0.5 day? 
Tmax Carrying capacity of T-cells 1500 mm~ 
N Bursting term for viral production after lysis 10-2500 
c Clearance rate of virus 2.4-3 day™! 
T Virus replication time 0-2 days 
Nrt Reverse transcriptase inhibitor efficacy (0,1) 
Np Protease inhibitor efficacy (0,1) 
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